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,^,: Abstract 

A new mechanism of longitudinal confinement of optical energy via coupled plasmon modes is proposed 
in chains of noble metal nanoparticles embedded in a graded dielectric medium, which is analogous to the 
^ ■ confinement of electrons in semiconductor quantum wells. In these systems, one can control the transmis- 
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sion of optical energy by varying the graded refractive index of the host medium or the separation between 
the nanoparticles to realize the photonic analogue of electronic transistors. Possible passband tunability 
by nanoparticle spacing and modulation of the refractive index in the host medium have been presented 



F^ ' explicitly and compared favorably with numerical calculations. 
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One of the most attractive studies in advanced optics is to conquer the general difficulties in con- 
ventional photonics: (1) restriction of the diffraction limit, which prevents achieving high degree 
of miniaturization and monolithic integration of optical devices and circuits; and (2) light fields 
existing in three dimensions, inhibiting integrated on-chip devices. A prominent way to overcome 
these two major obstacles is related to use of surface plasmons (SPs) in plasmonic structures, 
e.g., see Refs. 1-8 and references cited therein. Along with the advancement of nanofabrica- 
tion techniques, many nanostructures and devices based on SPs have been studied, for example, 
metallic nanoparticles,^ nanowires,- nanostrips,- and metal grooves^ and wedges,— as well as in- 
tegrated components such as nanoparticle chains,^ multilayers,— and combined plasmonic and 
dielectric waveguides. ^^ They essentially utilize SPs features of lateral/transverse localization and 
evanescence,'^ and their operations in near- field regime, thus can easily beat the above mentioned 
restrictions in conventional photonics. 

Although these plasmonic systems offer the advantages to relieve the difficulties in the minia- 
turization and integration of photonic devices, SPs suffer from intrinsic dissipation of consisted 
material, e.g., mostly noble metals like gold and silver and thus decay excessively along the prop- 
agation direction, or referred to the longitudinal direction. Nevertheless, their propagation length 
is still satisfactory in nanooptics, e.g., 1/e propagation lengths of a few tens of fiva are reported.— 
It is, alternatively, desirable to have localized SPs in the longitudinal direction, in addition to their 
naturally lateral/transverse confinement. In this sense, they may sustain modes in ultrasmall vol- 
ume and highly confined electromagnetic field. To this end, only very few works have focused 
on this issue, for instance, nano-concentration of optical energy by localizing SPs in graded plas- 
monic waveguides (i.e., axially nonuniform waveguides),— and symmetry breaking in plasmonic 
chains, which helps the enhancement of local fields, and facilitates the forming of localization of 
SPs due to asymmetric polarization charge distribution.— Other examples include using nanoscale 
Fabry-Perot cavity to confine SPs inside,— and constructing plasmonic crystals by decorating or 
structuring metallic surfaces at nanoscale. ^'^ In analogy of photonic crystals, coupled SPs show 
localization near the band-edges in these plasmonic crystals and demonstrate desirable localiza- 
tions. 

In this work, we show that precise control of longitudinal localization of SPs in plasmonic 
nanoparticle chains (PNCs) is possible by simply imposing a gradient in the refractive indices 
of the surrounding host medium, which can behave like a plasmonic switcher. The localization 
mechanism is different from those of the two extreme cases: periodically structured plasmonic 



crystal o^i^^i^" and disordered plasmonic systemsj^i'22, In fact, closely-separated PNCs were pro- 
posed for potential applications in energy guiding and routing inside the subwavelength scale. - 
Also they were demonstrated efficiently converting propagating light (far-fields) to SPs (near- 
fields).— The precise control of SPs in our proposed system can also achieve confinement of elec- 
tromagnetic fields in specific region to control their interactions with specimen.^"* Furthermore, we 
can envision applications of such PNC arrays that evanescently couple two dielectric slab waveg- 
uides, e.g., in photo-addressable birefringent polymer matrix (optical control) or in electro-optical 
materials such as liquid crystals (electrical control). 

Let us assume that a chain of N identical metallic nanoparticles of Drude form dielectric func- 
tion em{^) = 1 — ^p/^^ (lossless) is immersed in a dielectric host with its dielectric function 
varying from the left-hand side to the right-hand side along x-axis as 

e2{x) = eieft + ex /I, (1) 

where x denotes the position of the particles, / the total length of the chain, and c the coefficient of 
dielectric gradient which can be formed naturally, induced by a pump light that possesses intensity 
variation spatially, or by possible ultrasonic induced gradient in a piezoelectic host.— The i-th and 
the j-th nanoparticles are separated by a distance dij = d^li — j\, with their position denoted by 
Xi = do(i — 1), where i = 1,2,- ■ ■ , N. This means that the nanoparticles are equally spaced 
in a periodical chain which is immersed in a host with dielectric function increasing from eiett at 
the chain's leftmost site to eieft + c at the rightmost site. Therefore, each spherical nanoparticle 
acquires an isotropic dipole polarizability 

A = e.(x.)a3 ^-/"^7^("-\ . (2) 

em(uj) + 2e2 xj 



The dipolar resonant frequency is denoted by Ui = Juj^/[1 + 2e2(xj)] which makes the real part 
of denominator of Eq. Q zero.- In the coupled dipole approximation, the self-consistent coupling 
equations of the local fields Ej around the i-th nanoparticle read,- 

E, = Y,f{z,j)-i(3,E,) + Ef\ {t = l,2,---N), (3) 

where E^^ ' represents the external electric field, and the dipole coupling is given by T{i,j) = 
T{i,j)/e2{xi) where T{i,j) denotes the near field coupling in vacuum. In Cartesian coordinates 
(x, y, z), the (/3,7) components are given by 

rp (■ -N _ ^djj^pdjj^^ — \aij\ Ofj^y 
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Where djj represents the vectorial distance with /3 or 7 component dij^p^^. After substituting the 
Drude form into Eq. Q and eliminating the bulk plasmon frequency cjp in Eq. Q, we have the 
following coupled equation of motion for the dipole moments 
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where Al = 2 and Ax = — 1 are for the cases of longitudinally (i.e., parallel to the chain axis) 
and transversely polarized E^^ ' (simply indicated by E^ '), respectively. In view of the fact that 
we are concerned with the plasmon waves near the resonant frequency r^iSS i.e., cu ~ tUj, we have 
obtained Eq. dSj) with the approximation [e2(xj) — l]iu'^ + [2e2(x,;) + l]ujf ^ 3e2{xi)ujf. Equation 
^ captures the main features of the plasmonic coupling between nanoparticles in the PNCs. For 
example, it shows that the coupling depends on the resonant frequencies, the polarization, and the 
properties of the metal and the surrounding host, as well as the separations between the particles .- 
We seek plasmonic eigenmode solutions $q,(z) of Eq. ^ in the absence of the external driving 
term E^ ' . The results are shown in Fig. [T] for the longitudinal polarization and with parameters 
listed below: eieft = 3.0, c = 1.0, N = 100, and dg = 3a, namely / = 3a(N — 1), where a 
is the radii of the nanoparticles. In this nanoparticle separation, complications from multipolar 
interactions are negligible^^ and our approximations are quite accurate. It is noteworthy that in 
the calculations we used the approximation that the dielectric constant of the host around the i-th 
nanoparticle is homogeneous, e.g., denoted by e2(xj). In this case, we actually have a host dielec- 
tric function 62 (xj) gradually varying from 3.0 at the left-hand extremity to 4.0 at the right-hand 
extremity. The solid curve in Fig.lUa) represents the relationship (i.e., pseudo-dispersion relation) 
between the coupled plasmon frequency cu and the mode index a for the longitudinal polarization. 
Longitudinal results for the same PNC within homogenous hosts of dielectric constant 62 = 3.0 
(dash-dotted line) and €2 = 4.0 (dashed line) are also shown. It is noticed that the presence of 
a gradient in the host refractive index alters the plasmon (pseudo-) dispersion relation dramati- 
cally. This is also signified in the density of states (DOS) D(uj) = J2a ^(^ ~ ^a)/N, as shown in 
Fig.IHb). In detail, the DOS peaks at band-edges of the homogenous host cases (dashed and dash- 



dotted lines) are shifted toward the Mie dipole plasmon resonant frequencies ujq = ujp/^\ + 2e2, 
due to the gradient in the host refractive index. More specifically, for the graded host case we 
studied, there are two peaks in the DOS at around ujn = 0.331u;p and cufj = 0.364co'p, for near field 
full couplings (Full) between nanoparticles. The two peaks, however, appear at ujd = 0.337^;^ 
and ujc2 = 0.367u;p if only the nearest-neighboring couplings (NN) are taken into account. These 



two peaks in fact correspond to two points of transition between confined modes and extended 
modes, which is easily seen in Fig.lljc) showing the inverse participation ration (IPR) of the nor- 
mal mode,— as well as directly from the mode patterns (right panels in Fig.[l]). Interestingly, the 
modes outside the frequency range un < to < ujf2 are localized at the two extremities of the 
PNC, i.e., modes with to > uj{2 (high frequencies) are localized at the left-hand side while modes 
with to < ton (low frequencies) are confined at the right-hand side with dramatically increased 
IPRs. These localization/confinement are clearly shown by the mode patterns in Figs. [Hd) and 
1(f). However, normal modes with frequency ton < to < ^n are extended over the entire PNC 
[see Fig. [He)]. They exhibit IPRs that are near constant as shown in Fig.[llc). 

Intuitively, the existence of localized (confined) modes at the two ends is reasonable in view of 
the fact that the Mie resonances is lower at the right hand where there is a relatively larger refractive 
index in the host and higher at the left hand side with a relatively smaller refractive index in the 
host. But to further understand the localization-delocalization transitions, let us neglect long range 
interactions and truncate the summation in Eq. Q up to the nearest neighbors, i.e., consider NN 
couplings. The resulting coupled plasmon system is similar to a chain of graded coupled harmonic 
oscillators22 with additional on-site harmonic potentials. In this case, the particle mass nii, and the 
strength of the additional harmonic spring f/j are respectively 



1 + 2e2(x, 
3Xtofe2{xi 



m = IZTTi ^ {^ = 1,2,...,N), (6) 



U. = W?^?^ -2Ko, (. = 1,2,... ,iV), (7) 

3Ae2(Xi) 

where Kq = {a/dof" is a force constant between adjacent particles that depends on the ratio of the 
interparticle spacing and the particles radius. Figures |2t a) andEtb) show rrii and Ui for the case of 
cLq = 3a (solid lines) and do = 2.5a (dashed lines) respectively, for A = Al = 2. In this case rui 
and Ui are respectively increasing and decreasing functions of i, and for A = At = — 1, they are of 
the same functions of i (not shown). Illustrating this situation, we can understand that the lower 
frequency modes are confined at the right side of the chain and higher frequency modes at the left 
side. The lowest frequency cul (i.e., lower bound) of extended modes is the lowest frequency of 
the homogeneous system with m^ = mi and Ui = Ui for any i. This frequency wl is simply given 
by 



^L = VUi/m,. (8) 

Also, the highest frequency tou (i.e., upper bound) of extended modes is the highest frequency of 



the homogeneous system with rrij = rriN and f/j = Un for any i. This frequency is given by 



Wh='\/ • (9) 

V mN 

For parameters A^ = 100, cLq = 3a, c = 1.0, e^x = 3.0, and Al = 2, Eqs. (El) and ^ result 
in ciJl = 0.340u;p and cuh = 0.365ciJp, respectively. These are in excellent agreement with the 
numerical data of nearest-neighboring coupling as shown in Figs.IIJb) and 1(c), where cjl and cjh 
are indicated by the two vertical dashed lines. In Fig. Etb), they coincide respectively with the 
two peaks in the DOS with NN couplings, which locate right at the two transition frequencies of 
tUci = 0.337c<jp and oj^i = 0.367ujp. 

The analytic results of the transition frequencies given by Eqs. © and ® help to design the 
passband substantially, although they are for nearest-neighboring coupling. In this sense, we 
hereby attempt the tunable passband for the longitudinal polarization represented by the analytical 
forms and confirm it by numerical calculations. The results are demonstrated in Figs.|2tc) and 2(d), 
where we assume an additional nonlinear host dielectric variation, i.e., £2(0;) = (eieft + cx//)(l+p). 
Here p indicates the change of dielectric constant due to optical nonlinear effects, i.e., p = X-Epump' 
where x is the third order nonlinear susceptibility of the host media and -Epump represents a strong 
pump field. In Fig. |2tc) we plot cjl (solid lines) and cuu (dashed lines) as functions of c with 
p = 0.0 for do = 3a and do = 2.5a, respectively. We also show the results of u^ (symbol D) and 
ujc2 (symbol Q) extracted from the numerical calculations. The numerical data agree very well 
with the analytical ones. This figure indicates that with fixed nanoparticle spacing do in the PNC, 
cjl keeps relatively constant, while cjh decreases almost linearly for an increased gradient coeffi- 
cient c from, say, 0.2 to 1.5. On the other hand, increasing the separation of the nanoparticles, e.g., 
changing do = 2.5a to do = 3a narrows the passband dramatically, as both cjl and uu approach to 
the Mie resonant frequency cui because the particles become more isolated, that is, the couplings 
between them are diminished. 

Figure |2td) demonstrates tunability of the passband by nonlinearity effect in the host medium. 
For example we show the transition frequencies cjl and cuu versus p for do = 3a with two different 
gradient coefficient c = 0.5 and c = 1.0, respectively. Similarly the numerical data ofcud (symbol 
D) and coci (symbol Q) ^e compared favorably to the analytical results (curves). It is seen that for 
increased p, i.e., increased dielectric constant in the graded host, both the low and high transition 
frequencies are red-shifted. These are ascribed to the fact that the Mie resonant frequency is 
redshifted when the dielectric constant of the surrounding host is increased. Interestingly, u;l (cUci) 



decreases relatively at the same fashion as cjh (uJci) does, varying over a span, e.g., from around 
O.SbcUp to a extremely low frequency. Admittedly, it is hard to achieve such a large nonlinear 
modulation, for example, for p > 1.0. However, a relatively small modulation in the host dielectric 
nonlinearity, e.g., from p = 0.0 to p = 0.5 still results a desirable tuning of cjl and cjh- We believe 
these results can find potential applications of dynamical switching in plasmonics. The pump 
field may also induce a nonlinearity in the metallic particles. However, the nonlinearity in host 
is assumed to be larger, resulting in a graded polarizability along the chain. Finally, the results 
for the transverse modes (A = — 1) are quite similar to those of the longitudinal modes (A = 2) 
presented here. 

In summary, we propose the use of a gradient dielectric host to achieve the precise confinement 
of electromagnetic fields in plasmonic nanoparticle chains. It is demonstrated that the coupled 
surface plasmons in these systems become easily tunable between confined modes and extended 
modes. We offer an interpretation of the transition mechanism by mapping the problem onto a 
graded coupled harmonic oscillators. In addition to possible operation for plasmonic switching, 
these systems may be also of interest for applications such as surface-enhanced Raman scattering 
or biosensing. 
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Figure Captions 



FIG. 1: (Color online) (a) Pseudo-dispersion relations for full coupling. Those for cases of homogeneous 
host with €2 = 3.0 (dash-dotted line) and €2 = 4.0 (dashed line) can be regarded as the usual dispersion 
relations, (b) Density of states (DOS) versus coupled plasmon mode frequency for the full coupling. The 
line with circles represent results of graded case for nearest-neighboring coupling, (c) Inverse participation 
ratio (IPR) versus plasmon mode frequency, (d) Typical excitation in low frequency regime {uj = 0.315u;p). 
(e) Typical excitation near the resonant frequency {u: = 0.353a;p). (f) Typical excitation at high frequency 
{uj = 0.390a;p). 



FIG. 2: (Color online) (a) Effective mass rrii and (b) on-site potential Ui for do = 3a and cLq = 2.5a with 
c = 1.0 and p = 0.0. (c) Analytical (curves) and numerical (symbols) transition frequencies ljl (wh) and 
"^ci i'^ci) versus the host gradient coefficient c without nonlinearity (p = 0.0) for different nanoparticle 
distances d^. (d) Analytical (curves) and numerical (symbols) transition frequencies versus nonlinearity 
(pump field) modulation in the host medium, with d^ = 3a for different host gradient coefficient c = 0.5 
and c = 1.0, respectively. 



0.40 


-(a)' ' 
e =3.0 ^.- 


II 




' \v - 


^^""^^ 


0.35 




E, ^ 4.0 _ 




j/^ ^ 


2 


0.30 


^1 1 1 


c=1.0 _ 

1 1 



~^ I ^ \ ^ tT"" r~^ 
j;dj co=/0.315( 




8, =3 >8.. . =4~ 

- low I ^ highi 
J I I I I I I I l_ 



e 



-0.2 



1 20 40 60 80 1001 20 40 60 80 100 
Mode index a Site index i 



lo'f 



W) 



"m 




n — I — I — I — I — I — I — I — r 

Full ' j ^ U(e) CO = 0.353coJ o.2 




J I I I I I L 



e 



--0.2 



0.30 0.35 0.40 1 20 40 60 80 100 



co/co 



10 z 



1 _ \ CO 






10" 




«H ^ ^f) 



Site index i 

"1 — I — I — ] — I — I — r 




|C0 = 0.390co^ 0.2 



0.0 



J I L 



J I L 



e 



- -0.2 



0.30 0.35 0.40 1 20 40 60 80 100 



co/co 



Site index i 



Fig. 1 ./Xiao, Yakubo, and Yu 



10 




— do=3a (cj 



.®«, 



■GW, 



'Q>( 




d„ = 2.5a 



p = 



■"-^^ 



'-0, 



-©~ 



dddddddddddddd 



□DaDDDaDDD n aD n 



0.45 



0.40 3 



0.35 o^ 



0.30 



1 20 40 60 80 100 







0.25 S 



20 40 60 80 100 0.5 1.0 1.5 2.0 
i P 



Fig. 2. /Xiao, Yakubo, and Yu 



11 



